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The existence of an ^ = 1 avoided crossing in the spectrum of the solar-like pulsator CoRoT**-target HD 49385 was 
established by Deheuvels & Michel 2009. It is the first confirmed detection of such a phenomenon. The authors showed 
in a preliminary modeling of the star that it was in a post main sequence status. Being a 1.3A/Q-star, HD 49385 has had 
a convective core during its main sequence phase. The /i-gradient left by the withdrawal of this core bears information 
about the processes of transport at the boundary of the core. We here investigate the constraints that the observed avoided 
crossing brings on the /^-gradient in the core of the star. 
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1 Introduction 



2 Frequency at which the avoided crossing 
occurs 



Avoided crossings between non-radial oscillation modes have| 

been the object of many theoretical studies since |Scuflaire (197 jff evolution code CESAM^(^^^p_combined| 
found that low-order p and g modes could present a mixed 



with the oscillation code LOSC (IScuflaire et al. 2008] ). we 



behavior in condensed polytropes In particular, Dzi embowski d'm^mykW{\Wti^°'^^'^^ reproducing reasonably well 

both me classical and the mean seismic indi 



stressed that the frequency at which the avoided crossing 
occurs has a high diagnostic potential, especially to con- 
strain the processes of transport at the boundary of a convec- 
tive core. [Deheuvels & Michel (2009)| have shown that low- 
degree avoided crossings are associated with a distortion 
of the ridge of the corresponding degree. The authors ex- 
plained that this feature is the effect of a coupling between n 
modes during the phenomenon of avoided crossing and sug- 
gested that it could be used to contrain the /i-gradient in the 
inner parts of the star Such a distortion was observed in the 
solar-like pulsator HD 49385, based on the analysis of 137 
days of photometric data collected with the space telescope 
CoRoT (Deheuvels & Michel 2009, IDeheuveIs et al. 2010l l.l 
The authors found that it could only be explained by an 
I = \ avoided crossing and this enabled them to conclude 
that the star is in a post main sequence stage. 

We here study how the observation of this avoided cross- 
ing constrains the processes of transport in the center of 
HD 49385. We first focus on the frequency at which the 
avoided crossing occurs. We then show that the curvature of 
the I — \ ridge brings complementary information on the 
chemical composition at the edge of the core. 



indices of HD 49385.1 
We were then interested in selecting among these models, 
those which show an avoided crossing at the same frequency 
as the observed avoided crossing v'^^. To locate the avoided 
crossing, we found it convenient to use the frequency of 
the mixed mode which has an essentially g-mode behavior, 
since it varies much faster than the frequencies of the other 
modes. Based on Deheuvels et al. (2010), the fitted eigen- 
frequency of this mode for HD 49385 is = 748.60 ± 
0.23 /iHz. We therefore needed to ensure that v"^^ corre- 
sponds to the eigenfrequency of a g mode. 

One way of doing this would be to locate the avoided 
crossings in each of our models. The problem is that several 
modes are involved and they are not the same ones from 
one model to the other It is more convenient to link the 
frequency of the avoided crossing to the equilibrium quan- 
tities. In the frame of the asymptotic theory, the eigenfre- 
quency Vn^i of a g mode of high order n and low degree I is 
approximately given by the equation 

1/2 

dr (n-l/2)7r 
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where 
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L''=l{i+l), 

N is the Brunt-Vaisala frequency, and ri, r2 are the turning 
points of the mode in the g-mode cavity. The integral I{v) 
defined as 
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Fig. 1 Profile of the squared Brunt- Vaisala frequency in 
the inner region of a MS and a PoMS model of HD 49385 
(full line). The dashed line shows the contribution of the 
departure from adiabaticity and the dotted line, the contri- 
bution of the gradient of chemical composition. 
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Fig. 2 Variations in the integral liv"^^) during the evo- 
lution of a PoMS sequence of models for HD 49385 (thick 
solid line). The dashed lines indicate the the solutions of Eq. 
|4]for n — 1, 2, 3. The dotted hne represents the observed 
value of the mean large separation Ai/ = 56.3 /zHz. 



can therefore be used to determine whether or not v"^^ cor- 
responds to the eigenfrequency of a g mode. Indeed, it is the 
case when there exists an order n such that Eq.[T]is verified 



for u = ^ill hereafter use I{v2c) as a rough in- 

dicator of how the frequency of the avoided crossing com- 
pares in the model and in the observations. We checked a 
posteriori that the models which verify Eq. [T] indeed have 
an avoided crossing around 

2. 1 Evolution of / ( 1^°^^ ) with the age 

The Brunt- Vaisala frequency N can be written as a function 
of the following derivatives: 



dlnT 



'ad 



d\nP 



dlnT 
dlnP 



(91n/i 
51nP' 



(3) 



where P, T and fi correspond to the pressure, temperature 
and mean molecular weight, respectively. The subscript S 
indicates that the definition is valid for constant entropy. It 
is convenient to split the expression of A^-^ in two contribu- 
tions: 

^(Vad-V) + -^ 
tin tin 



(Vad - V) + -^V/i, (4) 

±p Lip 

where Hp is the pressure scale height and g the gravity in- 
side the star The left part of the second member in Eq.|4]cor- 
responds to the departure from adiabaticity. During the main 
sequence, this term is null in the convective core. When the 
star reaches the end of the main sequence, the nuclear re- 
actions stop in the core, which ceases to be convective. The 
temperature gradient becomes smaller than Vad and the left 
term increases. The right term in Eq.|4]is caused by the gra- 
dient of chemical composition which is left by the receding 
convective core. We show in Fig. [1] that both contributions 
are clearly identified both in the main sequence (MS) and in 
the post main sequence (PoMS) phase. 

As the star evolves past the end of the main sequence, 
the inner regions contract to trigger hydrogen burning in 



shell. Consequently the density increases in the core, which 
causes the factor gj Hp to also increase in Eq.H] It is there- 
fore obvious that the integral increases during the 
evolution of the star This can be seen for a PoMS model 
of HD 49385 in Fig. ID We note that, since the mean large 
separation /S.v monotonically decreases during the evolu- 
tion of the star, we used this quantity instead of the age as 
a parameter of evolution. This representation enables us to 
easily spot the models which reproduce both the observed 
value of the mean large separation and the frequency of the 
observed avoided crossing. 

2.2 Evolution of / ( v"^"" ) with the mass 

We now consider models for which the mean large separa- 
tion matches the observed one, and we study the influence 
of the mass. As can be seen in Fig. [3] the line of the TAMS in 
the HR diagram becomes closer to the line of iso- Aj/ when 
the mass increases. Therefore, for larger masses, the star 
is closer to the TAMS when it reaches the observed value 
of the large separation Aj^obs and the value of the integral 
I{v'^^) is expected to be smaller. We indeed observe for 
HD 49385 that the integral decreases with increasing mass 
(see Fig. [3]). With a given physics, we can find one mass for 
which both the value of the mean large separation and the 
frequency of the avoided crossing are reproduced. 

2.3 Evolution of I{v'^^) with the overshooting 

The models we here consider are in the PoMS stage. Core 
overshooting is not supposed to play a direct role since the 
convective core has vanished. However, the profile of the 
gradient of chemical composition which was left behind 
as the convective core receded during the main sequence 
greatly depends on the amount of overshooting which ex- 
isted at the edge of the core. As is usually done, we im- 
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plemented overshooting in our models as an extension of 
the convective core over a fraction aov of the pressure scale 
height Hp. We studied the influence of aov on the frequency 
of the avoided crossing. 

There are two contradictory effects of core overshooting 
on the integral 

- First, by bringing more hydrogen to the core, the over- 
shooting enables the star to stay longer in the MS stage. 
This can be seen in Fig. [3] in the case of HD 49385. 
Consequently, for a given mass, the bigger ctov, flie less 
the star has evolved from the TAMS when it reaches 
the observed value of the mean large separation. As ex- 
plained above, we therefore expect /(^'°^^) to decrease 
when Q!ov increases. We note that above a certain value 
of aov, the model already has a mean large separation 
smaller than Ai^obs when it reaches the TAMS, and no 
satisfactory model can be found. This can be used to 
set an upper limit on aov With the physics we use, we 
found that aov < 0.16 for HD 49385. 

- Secondly, when adding overshooting, the mixed core 
during the MS phase is bigger and so the resulting He 
core at the TAMS is bigger. Therefore, the inner parts 
of the star need to be denser to start burning hydro- 
gen in shell. In this case, the factor g/Hp in Eq.|4] in- 
creases, and the integral increases. This implies 
that should increase with aov 

For a mild overshooting, the star is far from the TAMS when 
it reaches Ai/obs, and the first effect has little influence. We 
can see in Fig. [3] that I{v2c^) increases with aov for low 
values of aov. This tendency is reversed for large values 
of aov Indeed, with a larger amount of overshooting, the 
star is very close to the TAMS and the first effect prevails. 
Due to the competition between these two effects, certain 
masses lead to two different models fitting both Ai^obs and 
the frequency of the avoided crossing: one with a mild over- 
shooting away from the TAMS, and one with a larger aov, 
close to the TAMS. This is the case for M — 1.36M0 (see 
Fig.H. 

In Fig. m all the models which lie on the dashed line 
fit the position of HD 49385 in the HR diagram, the ob- 
served value of the mean large separation, and the frequency 
at which the avoided crossing occurs. We found that aov G 
[0, 0.16] but at this stage the overshooting parameter cannot 
be further constrained. 
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Fig. 3 Top: Evolutionary tracks in the HR diagram (left) 
and values of the integral for models of different 

masses fitting At^obs for HD 49385, with given physics 
(aov = 0). The masses range from 1.3Mq (light gray) to 
I.SGA/q (dark gray). In the HR diagram, the box represents 
the position of the star within l-cr error bars and on the dot- 
ted line, Ai' — Aj/obs- In the right plot, the dashed line 
corresponds to the solution of Eq. 2] for n = 2. Bottom: 
Same as the top, but for models with different aov and fixed 
mass {M = I.3IM0). The overshooting parameter ranges 
from (light gray) to 0.15 (dark gray). 



4.5 




3.0 - 

3 Curvature of the ridge 



For low-degree modes, the p-mode cavity and the g-mode q qq q q ^ q .| q q 

cavity are strongly coupled since the evanescent zone that a„ 
separates them is small. [Peheuvels & Michel (2009)| showed[ 

that in this case, avoided crossings involve more than two Fig. 4 Variations in the integral with overshooting 
modes: for one mode trapped mainly in the core, several for different masses. The masses range from 1.3Mq (light 
consecutive-order modes trapped mainly in the envelope have| gray) to l.SOA/© (dark gray). The dashed fine corresponds 
a non-negligible part of their kinetic energy in the g-mode to the solution of Eq.|4]for n = 2. 
cavity. Their frequencies are modified, resulting in a charac- 
teristic distortion of the ridge, whose shape depends on the 
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size of the evanescent zone. We used the curvature of the 
ridge as an additional constraint for our models of HD 49385.| 

3.1 Comparison with the observations 

To efficiently compare the curvatures of the ridges, the fre- 
quency at which the avoided crossing occurs in the model 
must match precisely that of the observations. The best way 
to ensure this is to check that the mode which is trapped 
mainly in the core has the same frequency in both spectra. 
Indeed, its frequency varies much more rapidly than that of 
the other modes. Until now, we have used /(t'"^^) to locate 
the avoided crossing. As stated above, it is only a crude in- 
dicator. Here, we made small adjustments in the age of the 
models selected in Sect. [2] so the mode which has a dom- 
inant g-mode behavior has the frequency of the observed 
mode (mode tti in Deheuvels et al. 2010). We then used the 
individual £ = 1 large separations to compare the curvatures 
of the ridges. We found that for all the satisfactory models 
of Sect. |2] the curvature of the £ = 1 ridge is less pro- 
nounced than in the observations (see Fig. |5]l. This means 
that in these models, the coupling between the g-mode cav- 
ity and the p-mode cavity is underestimated. It was pointed 
by [Shibahashi & Osaki (1976)] that the strength of the cou- 
pling increases when the size of the evanescent zone de- 
creases. The size of this region needs to be decreased in the 
models. 

The outer turning point of the g-mode cavity and the 
inner turning point of the p-mode cavity constitute the lim- 
its of the evanescent zone. While the latter varies very little 
in the different models we considered, the former strongly 
depends on the /^-gradient which was left behind by the re- 
ceding core during the main sequence. Indeed, we can see 
in Fig. [U that this point is located in the V/i component of 
the Brunt- Vaisala frequency. 

3.2 Effect of overshooting 

Adding core overshooting has the effect of extending the 
size of the mixed region associated with the convective core 
and the /x-gradient left by the withdrawal of the core is lo- 
cated less deeply in the star. We could therefore expect over- 
shooting to increase the coupling, as required. However, we 
observe that for HD 49385, the evanescent zone is even 
wider with overshooting and the agreement with the obser- 
vations for the curvature of the £ ^ I ridge is worse (see 
Fig. |5]l. The case of HD 49385 is particular since the star 
is close to the TAMS. The models without overshooting are 
evolved enough after the TAMS so that the nuclear reac- 
tions in shell have smoothed the /x-gradient. This has the 
effect of reducing the size of the evanescent zone. On the 
contrary, the models with overshooting reach Aj^obs very 
close to the TAMS and the /i-gradient is exactly the one 
which was left when the core receded. And so, eventhough 
the /i-gradient is located less deeply, the evanescent zone 
is not smaller This raises the question of the treatment of 
overshooting. Here, in the overshooted region, the mixing is 
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Fig. 5 Left: Propagation diagram of two models of 
HD 49385 fitting the position of the star in the HR dia- 
gram, the observed Av and the frequency of the avoided 
crossing: one without overshooting (full lines) and one with 
Oov — 0.7 (dashed lines). The Brunt-Vaisala frequency is 
shown in black and the Lamb frequency for £ ^ 1 modes 
in gray. The dotted line indicates ly"^^. Right: Profile of the 
£ = 1 large separation for the same models. The filled cir- 
cles correspond to the obervational values with l-a error- 
bars. 

considered total and the gradient is supposed to be adiabatic 
(known as instantaneous mixing). Such a simple treatment 
of overshooting is found incompatible with the observations 
forHD 49385. 

4 Conclusion and perspectives 

Two features of the £ = 1 avoided crossing observed in 
HD 49385 were reckoned to give constraints on the /i-gradient§ 
in the inner parts of the star: the frequency at which the 
avoided crossing occurs and the curvature of the £ — 1 
ridge around the avoided crossing. By computing evolution 
models of HD 49385 fitting the global parameters, the mean 
seismic indices and the frequency of the avoided crossing, 
we showed that only an upper limit to the amount of over- 
shooting that existed at the edge of the core can be set: 
Qfov < 0.16. Among all these models, none is able to sat- 
isfactorily reproduce the strong curvature of the observed 
£ ^ 1 ridge. We conclude that the crude description of over- 
shooting as an instantaneous mixing which we used here 
(and which is most commonly used) is incompatible with 
the observed spectrum of HD 49385. We plan to test other 
descriptions of the mixing which exists at the edge of the 
convective core. For instance, this mixing could be treated 
as a diffusive process jVentura et al. 19981 ). 

Acknowledgements. 



© 2006 WILEY- VCH Verlag GmbH & Co. KGaA, Weinheim 



www.an-joumaLor^ 



Astron. Nachr. / AN (2006) 



5 



References 

Deheuvels, S., & Michel, E. 2009, Ap&SS , 241 

Deheuvels, S., et al. 2010, arXiv: 1003.4368 

di Mauro, M. P., Christensen-Dalsgaard, J., Patemo, L., & 

D'Antona, F. 2004, Sol. Phys., 220, 185 
Dziembowski, W. A., & Pamyatnykh, A. A. 1991, Astron. Astro- 

phys., 248, LU 

Mathis, S., Palacios, A., & Zahn, J.-P. 2004, Astron. Astrophys., 
425, 243 

Miglio, A., Montalban, J., & Maceroni, C. 2007, MNRAS , 377, 
373 

Morel, P. 1997, Astron. Astrophys. Suppl. Ser., 124, 597 

Osaki, J. 1975, PASJ , 27, 237 

Scuflaire, R. 1974, Astron. Astrophys., 36, 107 

Scuflaire, R., Montalban, J., Theado, S., Bourge, P.-O., Migho, A., 
Godart, M., Thoul, A., & Noels, A. 2008, Ap&SS , 316, 149 

Shibahashi, H., & Osaki, Y. 1976, PASJ , 28, 199 

Ventura, P., Zeppieri, A., Mazzitelli, I., & D'Antona, F. 1998, As- 
tron. Astrophys., 334, 953 



www.an-joumal.org 



© 2006 WILEY-VCH Verlag GmbH & Co. KGaA, Weinheim 



